Imperial/TP/03-4/7 

Disentanglement and Decoherence by Open System Dynamics 

P.J.Dodd and J.J.Halliwell 

Blackett Laboratory 
Imperial College 
London SW7 2BZ 
UK 

Abstract 

The destruction of quantum interference, decoherence, and the destruction of entanglement 
both appear to occur under the same circumstances. To address the connection between these 
two phenomena, we consider the evolution of arbitrary initial states of a two-particle system under 
open system dynamics described by a class of master equations which produce decoherence of each 
particle. We show that all initial states become separable after a finite time, and we produce the 
explicit form of the separated state. The result extends and amplifies an earlier result of Diosi. We 
illustrate the general result by considering the case in which the initial state is an EPR state (in 
which both the positions and momenta of a particle pair are perfectly correlated). This example 
clearly illustrates how the spreading out in phase space produced by the environment leads to 
certain disentanglement conditions becoming satisfied. 

PACS numbers: 03.65.-w, 03.65.Yz, 03.65.Ud 
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I. INTRODUCTION 



Much effort has recently been devoted to understanding the properties of entangled quan- 
tum states. This effort is largely driven by the emerging field of quantum computation and 
in particular the desire to manipulate entangled states in a practically useful way. How- 
ever, another reason why the study of entanglement is of interest concerns the question of 
emergent classicality from quantum theory. Entanglement represents the possibility of cor- 
relations which are greater than those anticipated in classical theories. Hence, any account 
of emergent classicality must explain how entanglement is lost. The explanation of this is 
in fact reasonably simple and is closely related to decoherence, the destruction of interfer- 
ence. The purpose of this paper is to discuss the destruction of entanglement in some simple 
systems and its connection to decoherence. 

A state of a bipartite system is said to be separable (or disentangled) if it may be written 
in the form 

where Pi > [lj. Such a state describes essentially classical correlations and can never violate 
Bell's inequalities. However, it turns out to be surprisingly difficult to determine, in general, 
whether a state may be written in this form. Peres [2j made the very useful observation 
that a separable p remains a density operator under the operation of partial transpose 
(transposition of one subsystem only). Hence, a necessary condition for separability is that 
density operator properties are preserved under partial transpose. This condition was shown 
by Horodecki to be sufficient in the case of 2 x 2 and 3x3 dimensions |3j], but generally not 
otherwise. In the case of continuous variables the Peres- Horodecki condition has a useful 
expression in terms of Wigner functions, where becomes 

W{p 1 ,x 1 ,p 2 ,x 2 ) = Pi W l A (p 1 ,x 1 )W l B (p 2 ,x 2 ) (2) 

i 

The Peres-Horodecki condition is then that W(pi, x±,p 2 , x 2 ) remains a Wigner function 
under p 2 — > —p 2 . It has been shown by Simon that the condition is both necessary and 
sufficient for the case when the Wigner function of a bipartite system is Gaussian . Duan 
et al. considered a different necessary and sufficient condition for the separability of 
bipartite Gaussian states based on the variances of a class of pairs of commuting EPR-like 
operators (of which X\ — x 2 and p\ +p 2 are an example). (See Ref. p for a discussion of the 



connection between these condtions). There are undoubtedly more conditions for Gaussian 
states. 

A closely related idea is that of entanglement-breaking maps. This is a map $ acting on 
a subsystem A such that (1^ ® $)(r) is separable for all choices of state r on A ® B, and 
for all finite-dimensional choices of B. A theorem due to Horodecki and Shor then states 
that a map $ is entanglement-breaking if and only if it has the Holevo form 

Hp) = ^2 R k Tr(F k p) (3) 

k 

where Fu are a set of POVMs and are are set of density operators which are independent of 
p 0,0,0 • Note that this result refers to the dynamics of one of the subsystems only. What is 
particularly interesting about this type of map, is that appears naturally in the open system 
master equations of the type frequently used in decoherence studies. In particular, Diosi 
[lol | has recently considered the open system dynamics described by the master equation 

i D 

p= -[H,p]-—[x,[x,p]] (4) 

This equation describes a particle coupled to a heat bath in the limit of high temperature 
and negligible dissipation and is the simplest equation used to describe decoherence. If we 
write the solution to this equation as 

Pt = Hpo) (5) 

then Diosi has shown that this map becomes entanglement-breaking after sufficient time has 
passed for the P-function to positive 1Q |. The time taken for this to happen is typically 
very short and is essentially the same as the timescale required for decoherence. Although 
what is particularly interesting is that complete disentanglement occurs after a finite time, 
unlike decoherence which, in the usual view of it, is asymptotic in time. (See also the similar 



result by Diosi for spin systems in Ref. 
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It is not hard to see that both decoherence and disentanglement tend to be produced 
under the same circumstances. A simple illustration of that fact is as follows. Consider first 
a one-particle system in an initial superposition state 

|V> = |V) + 10) (6) 

and thus with density operator 

p= 1^)^1 + 1^(01 + 10)^1 + 10)^1 ( 7 ) 



Suppose now it is subject to evolution according to the master equation Eq.(J3J), with solution 
written in the form Eq.fjSj). It is generally known that, if the initial state © is a superposition 
of localized position states, evolution according to (jlj) tends to kill the off-diagonal terms. 
That is, we have 

<WX0l) « (8) 

after a typically very short time. This means that the density operator becomes essentially 
indistinguishable from the evolution of the mixed initial state, 

pl=W)W\ + |0X0| (9) 

This is the simplest account of decoherence of a single particle coupled to an environment. 
Now suppose we consider a two-particle system in the entangled state, 

m = iv>i> ® iv>2> + \<h) ® I0 2 ) (io) 

with density operator 

p = IV'iXV'il ® h^X^I + l^i)(0il ® l^ 2 )(0 2 | 

+ |0l)(^l|® |02}(^2| + |0l}(0l|® \<h)(<h\ (11) 

If we now let both particles evolve according to the dynamics $ ® we find that once again 
the off-diagonal terms go away, so for example, 

<WiX0i|) « (12) 
and the density operator becomes indistinguishable from that obtained by the initial state, 

P = MM ® |^XV>2| + I0i)(0i| ® |0 2 )(0 2 | (13) 

which is separable. Hence, the mechanism that destroys interference also destroys entangle- 
ment. 

Another way to see why coupling to an environment will destroy entan glem ent is to 
appeal to the fact that the property of entanglement has an exclusive quality [12J. Suppose 
A is entangled with B. Then if B becomes entangled with a third party C it diminishes its 
entanglement with A. Hence an environment coupling to one or both of the two particles in 
an entangled state will cause one or both of them to become entangled with the environment, 
thereby diminishing their entanglement with each other. 

4 



The aim of this paper is to investigate the destruction of entanglement through interacting 



with an environment, extending and elaborating the earlier result of Diosi 
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lll | and others 



14] 



In Section 2 we consider the dynamics of a particle coupled to an environment, concen- 
trating on dynamics of the form Eq.(j3J). This is reasonably standard material but we write 
it in a form which is most useful for studying disentanglement. In Section 3, we consider the 
evolution of bipartite systems under the dynamics of Section 2. We show that an arbitrary 
initial state achieves the explicitly separated form (J2J) after finite time. In Section 4, we 
consider the evolution of the EPR state in the presence of an environment. This simple 
example gives a clear picture of how the various separability conditions come to be satisfied 
as a result of interacting with the environment. We summarize and conclude in Section 5. 



II. EVOLUTION IN THE PRESENCE OF AN ENVIRONMENT 

Before considering the evolution of entangled states in the presence of a thermal environ- 
ment, it is useful to consider first the simplest case of a single particle coupled to a thermal 
environment in the limit of high temperature and negligible dissipation, with no external 
potential. The master equation (JU) for the density matrix p(x, y) is, 

dp ih ( d 2 p d 2 p\ D 



dt 2m\dx 2 dy 2 J h 2 ^ ^ 9 (M) 
where D = Im^kT . In the Wigner representation, the corresponding Wigner function 

W(p, x) = ^ J dt e-** p(x + \(L,x- (15) 

obeys the equation, 

dW _ ^p_dW_ D ^W_ 
Following Ref.|15|. this equation may be solved in the form 

W t (p, x) = J dpodxo K(p, x, t\p , x , 0) W (p , x ) (17) 
where K(p, x, t\p , x , 0) is the Wigner function propagator, and is given by 

K = exp [-a(p - p d ) - (3(x - x d ) - e(p - p d )(x - x d )] (18) 
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where p c i and x c i denote the classical evolution from p , xo to time t, 

Pel = PO, Xcl = x -\ (19) 

m 

(For convenience, we ignore exponential prefactors unless necessary). The coefficients a, (3 
and e are given by 

1 „ 3m 2 3m , , 

In fact, the general form of the propagator (|T5j) can be used to describe the most general 
type of linear dynamics - arbitrary environment temperatures, non-negligible dissipation 
and the inclusion of a harmonic oscillator potential - for suitable choices of a, /3, e and p c i, 
%d Ql- More general dynamics are considered in another paper 

With the simple change of variables Xq — > x$ — Pot/m we may write, 

W t (p,x) = J dpodxo exp[-a(p-po) —/3(x — x ) — e(p — p )(x — x )] 

x W^po,x ) (21) 

where 

Wq(p 0} x ) = W (p 0l x - pot/m) (22) 

This simple transformation is a linear canonical transformation, which corresponds to a 
unitary transformation of the initial state, so Wq is still a Wigner function. For decoherence 
and disentanglement, the important aspects of the evolution are contained in the convolution 
with the exponential function. Following Diosi and Kiefer [l?j], it is now very useful to 
introduce the notation 



z 



P 

y x 

and also to introduce a class of Gaussian phase space functions 

1 



(23) 



g(z;C) = exp^--z T C- 1 zj (24) 

The 2x2 matrix C is positive definite and \C\ denotes its determinant. The phase space 
function g(z\ C) is a Wigner function if and only if 

\C\ > \ (25) 

(This is essentially the uncertainty principle). A useful result is the simple convolution 
property, 

J d 2 z g{z l - z; C)g{z - z 2 ; B) = g{z l - z 2 , C + B) (26) 



We can use these Gaussians to compute smeared Wigner functions by convolution: 

W(z) = J d 2 z' g(z - z'; C)W{z') (27) 

Then it follows that the smeared Wigner function will be positive if and only if (|25jl holds. 
This is because the smeared Wigner function is then equal to the overlap of two Wigner 
functions, for which we have the result 

J dpdx W P1 {p, x)W P2 (p, x) = — Tr ( PlP2 ) (28) 

which is clearly always positive. For example, the Q-function, which is always positive, is 
obtained in this way by smearing with a minimum uncertainty Wigner function. In this 
notation the propagation Eq.(|21|l of the Wigner function is 

W t (z) = J d 2 z' g{z - z'; A) (29) 

For the free particle without dissipation considered above we have 

/ 2 t/m \ 
A = Dt\ (30) 
\t/m 2t 2 /3m 2 l 

and therefore 

D 2 t A 

showing that the Wigner function tends to spread out with time. 

Using the above description of the dynamics, it is straightforward to compute the vari- 
ances of x and p after a time t. They are 

(Ap) 2 = 2Dt+(Ap) 2 (32) 

(Ax), 2 = ~ + (Ap)l— 2 + -a(x,p) + (Ax) 2 (33) 
3 m z m z m 



where 



evaluated in the initial state. 



P) = j ( £ P + px) - (x) (p) (34) 
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III. EVOLUTION OF BIPARTITE STATES IN THE PRESENCE OF AN ENVI- 
RONMENT 



Consider now the case of a two-particle system in an initially entangled state. The 
two particles are not coupled to each other, but are each separately coupled to a thermal 
environment, as described in the previous Section. 

For our two-particle system, the Wigner evolution equation for the Wigner function 

W(p 1 ,x 1 ,p 2 ,x 2 ) = W(z 1 ,z 2 ) is, 

dW _p ± d^_p 2 dW_ D &W_ D d 2 W 
dt m dx\ m dx 2 dp 2 dp 2 

This equation may again be solved using propagators (with the unitary part removed, as in 
Eqs. (J2H), (1221) )• The two-particle Wigner function at time t is then given by 

W t (z u z 2 ) = J d 2 z[d 2 z' 2 g(z 1 - zi; A)g(z 2 - z' 2 ; A) W^z[, z' 2 ) (36) 

where the matrix A is given by Eq. (j30|) . We will show that this evolves into the explicitly 
separable form (J2J) after sufficient time has elapsed. The key idea is to write the matrix A 
in the propagator as 

A = Ai/4 + B (37) 

where -A1/4 is the matrix of variances for a minimum uncertainty Wigner function (whose 
explicit form will be given below). We use the convolution property to write, 

g(z x - zi; A) = J d 2 z x g(z x - z x \ B)g(z 1 - z[; A 1/4 ) (38) 

and similarly with g(z 2 — z' 2 ; A). On the right-hand side, the second Gaussian is a Wigner 
function because \Ai/±\ = h 2 /4. The first Gaussian with variance matrix B will also be a 
Wigner function if and only if 

\B\ = \A-A xfA \>^ (39) 
We may now write the two-particle Wigner function at time t as 

W t (z 1 ,z 2 ) = J d 2 Zl d 2 z 2 g{z 1 -z 1 ;B)g(z 2 -z 2 ;B)W {z 1 ,z 2 ) (40) 



where 



Wo(zi, z 2 ) = J d 2 z[d 2 z' 2 g{z x - z[; A 1/A )g{z 2 - z' 2 ; A 1/A ) W^z[, z' 2 ) (41) 



and is positive, as explained above. This is the desired result. The Wigner function at time 
t is of the separable form (J2J) as long as the matrix B satisfies (J3"9*j) . This is because the two 
Gaussians in (J4U)) are Wigner functions if Eq. (|39|) holds, and the smeared Wigner function 
Wo is positive and corresponds to the term pj in Eq. (]%]). 

We may also compare with the closely related result of Diosi [10], who showed that a 
one-particle system achieves the entanglement-breaking form (J3J) after finite time, and hence 
cannot be entangled with anything else. Tracing over particle 2 in Eq. (j40j) . we obtain 

W t (zx) = J d 2 Zl g( Zl - z i; B) Wo(zi) (42) 

where Wq is the smeared Wigner function of the one-particle system. This is clearly the 
Wigner transform of an expression of the form Q because firstly, the Gaussian g(zi — zj.; B) 
is a class of Wigner functions independent of the initial state so corresponds to Rk, and 
secondly, Wo is a Q-function so is of the desired form Tr (Fkp) with Fk taken to be a coherent 
state projector. Hence we completely agree with Diosi's result. However, this result takes 
Diosi's result further in that it gives the explicitly separated form of the disentangled state. 

Consider now the timescale on which the disentanglement condition Eq. (j39|) becomes 
satisfied, where the matrix A is given by Eq. (JHOJ) and Ax/4 is stil 
of this timescale was reported (without explicit details) in Refs. 
this analysis in a small way. We also need to compare the calculation with a result of the 
next section, so we give some of the details. 

Since the calculation involves a small numerical calculation, it is useful to define dimen- 
sionless variables i, p, x, defined by 

*=C^V*- (43) 



to be chosen. Calculation 
but we can extend 



D J 

p = (hmD)^ p (44) 

x = UcJ * (45) 

In these dimensionless variables, we now take the matrix Ax/4 to be 

A = H /V ~ 2 1/2 ) (46) 
^ 1/2 a^/V2j 

Refs [l7 1 made the particular choice <Jq = \pl on the grounds that this give a particularly 
robust evolution for certain sets of initial states [18]. For the moment, we leave it general. 



The condition Eq.()39|) now becomes, 

It3_2^2 f _l >Q (47) 

3 3 s ~ K ' 

where s = o\j\[2. The time at which this condition becomes satisfied may be estimated 
by plotting the function in Maple for various values of s. For the choice a$ = \/2 we have 
s — 1, and it is straightforward to then show that the the condition is satisfied for values of 
t greater than about 1.39. This means, 

*>1.97xC^V (48) 



2D J 



in agreement with Refs. 
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17| . (The factor of 2 is included because Refs.|l0l. Il7| used D/2 
to denote what we denote by D). 

However, if we tune the value of <jq to make the timescale as short as possible, then 
numerical experiments show that the optimal value is about s = 0.9 which gives the slightly 
shorter timescale 

t >1.9Bx(§* (49) 

This improvement is clearly insignificant, but it does however show that the Diosi-Kiefer 
choice is very nearly optimal. 

IV. THE EPR STATE 

Although the above results show disentanglement for general initial states, it is of interest 
to look in detail at a particular entangled state to see exactly how the entanglement goes 
away. We therefore consider the entangled state first introduced by Einstein, Podolsky and 
Rosen Instead of the two canonical pairs pi,x± and P2, £2 it is very useful to use instead 
the rotated coordinates, 

X = xi-x 2 , K = -{p 1 -p 2 ) (50) 
Q = ( Xl + x 2 ), P = ^( Pl +p 2 ) (51) 

This is a canonical transformation, with the new canonical pairs being K,X and P,Q. 
However, we also have the important relations, 

[X,P] = 0, [Q,K]=0 (52) 
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This simple observation is the basis of the EPR state, since it means we may choose a 
representation in which X and P are definite. In fact, EPR defined the state 



V(X,P) = 5(X)5{P) (53) 

This state is not strictly normalizable. We therefore instead consider Gaussian states which 
may be made arbitrarily close to this state. In particular, we consider a Gaussian state, 
which, in the rotated coordinates has Wigner function 

( K 2 X 2 P 2 O 2 \ 
W(K,X,P,Q)^^--.-^--.-^ (54, 

The widths ax and op may be chosen to be arbitrarily small, but for this to be a Wigner 
function, the remaining widths must satisfy 

h 2 h 2 

o K a x > — , a P a Q > — (55) 

(Because the transformation from the original coordinates to rotated ones is a linear canoni- 
cal transformation, it preserves Wigner function properties, so the conditions to be a Wigner 
function have the same form as in the original coordinates.) 

As an aside, we remark that Bell has noted that Gaussian states have a positive Wigner 
function, and therefore have hidden variable interpretation, so cannot violate Bell's in- 
equalities j^J. This perhaps suggests that they do not have any significant entanglement 
properties. This is, however, in terms of observables local in phase space quantities. It has 
been pointed out that there are other observables, in particular the parity operator, in terms 
of which Gaussian states do violate Bell's inequalities [2]]]. Hence it is of interest to study 
entanglement of Gaussians. 

In te rmS of the variables «, Q,, the Peres-Horodedd QBQ condition f 0r disentan- 
glement is the condition that W(K, X, P, Q) must remain a Wigner function when P and 
K are interchanged, that is, that 

( P 2 X 2 K 2 O 2 \ 
W(K,X,P,Q) = ^-—- — - — -^ (56) 

is a Wigner function. The conditions for this to be a Wigner function are, 

h 2 ft 2 
(r x °p > -J, °q°k > -J (57) 
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which is different to (|55|). In particular, the first of these relations will not be satisfied 
when ax and o> are chosen to be very small, which is the EPR case. This is generally 
the case for entangled states - they fail to remain Wigner functions under interchange of P 
and K because they are then too strongly peaked in phase space an violate the uncertainty 
principle. 

To watch the destruction of entanglement, we will work with the disentanglement condi- 
tion put forward by Duan et al. rather than the Peres- Horodecki condition. They wrote 
down a class of necessary and sufficient conditions for a Gaussian state to be disentangled. 
In terms of dimensionless variables K,X,P,Q (defined as in Eqs. (J3E|)-(j2H!0> one of those 
conditions is 

(AX) 2 + 4(AP) 2 > 2 (58) 

This is clearly not satisfied for the EPR initial state Eq. (J5%J| , but it is of interest to see how 
it becomes satisfied under evolution in the presence of an environment. The key point here 
is that the condition Eq. (j58j) means that the state is reasonably spread out in phase space, 
and this phase space spreading is precisely what evolution in the presence of an environment 
produces, so we expect that the condition will become satisfied after a short period of time. 

In the presence of an environment, the initial state will evolve according to Eq.(|35|). 
However, since the regularized EPR state (J53j) factors in terms of the rotated coordinates 
(J50|) . (fo"Tj) . it is more useful to use those, and in terms of them, the Wigner equation is 

~dT ~ ~dQ ~ ~dX + ~d? Y+ ~dK* ^ ' 

Interestingly, the dynamics also factors in the rotated variables and is essentially the same 
as the single-particle dynamics. Using Eqs. (j32"|) and (J33|) . it is then easily seen that 



(AP) 2 = 2Dt + a 2 P (60) 
8Dt 3 4 
3 m 2 m 2 



(AX) 2 = + ^ + 4 (61) 



Inserting these expressions in the condition Eq. (j58|) . we find that, in dimensionless variables, 
it is satisfied as long as 

8 1 

-f 3 + 4c? + 8t + — > 2 (62) 
3 4c 

Here, c is the dimensionless form of a\ and we have used the uncertainty principle to 
eliminate a\. The interesting case is that in which <j\ is very small, but we cannot set it 
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to zero because then <j\ would be infinite. We have also set crf> = 0, which represents the 
most extreme case. (We can do this because <Tg, its conjugate width, does not appear). 

The polynomial in Eq.(|62j) may be plotted using Maple for various values of c. We have 
found from these plots that for all values of c the condition becomes satisfied for values of t 
greater than about 0.19, that is, for 

i 

t -°- 27x i^y (63) 

This is the time after which the EPR state becomes disentangled, according the condition 
Eq.(|58|). This time is actually a lot shorter than the timescale Eq. ()49|h hence is safely 
consistent with the general result. The timescale Eq. (j49jl is the time for any initial state to 
become disentangled whereas the timescale Eq. (J63j) is only for the EPR state. This actually 
indicates that the EPR state is perhaps not a very entangled state, at least in terms of phase 
space variables. 



V. DISCUSSION 



We have shown, in a variety of different ways, how open system dynamics destroys en- 
tanglement in a system consisting of two entangled particles. Entanglement is destroyed by 
the same mechanism that destroys interference. In particular, we have shown that, under 
a simple open system dynamics, any initial two-particle state achieves the explicitly disen- 
tangled form Eq.(j2J after a short, finite time. We illustrated this general result with the 
particular case of the EPR state. 

The dynamics employed here are those of the free particle coupled to a bath of oscillators, 
in the limit of negligible dissipation. It is of interest to extend the analysis to include 
dissipation, a potential, and also to include interactions between the particles which will then 
compete with the disentangling effects of the environment. This is considered in another 

3- 



paper 



As previously noted 
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ll(, it is striking that complete disentanglement is achieved in 



finite time, whereas the interference terms in the density matrix tend to take infinite time 
to completely go away. That is, one is inclined to say that complete decoherence takes an 
infinite time. Although note that the word "decoherence" is used in a variety of different 
ways. 
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This perhaps suggests that the exact vanishing of the off-diagonal terms in the density 
matrix is perhaps too strong a condition for determining when a quantum system is essen- 
tially classical. It is, for example, often suggested that positivity of the Wigner function 
is a useful condition characterizing quasiclassicality, because this then means that the one- 
particle system is a hidden variables theory. Significantly, the Wigner function typically 
becomes positive after a finite time in open system dynamics. Hence, a reasonable meaning 
to attach to the word decoherence is that it is the situation under which prediction for all 
variables may be described by a hidden variables theory. These ideas will be explored in 
future publications. See Ref.|22j for related discussions. 
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